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The formation of breakdown pattern on an insulating surface under the influence of a transverse magnetic 
field is theoretically investigated. We have generalized the Dielectric Breakdown Model (DBM) for the case 
of external magnetic field. Concept of the Magnetic Fractal Dimensionality (MFD) is introduced and its 
universality is demonstrated. It is shown that MFD saturates with magnetic fields. The magnetic field 
dependence of the streamer curvature is obtained. It is conjectured that nonlinear current interaction is 
responsible for the experimentally observed 'spider-legs' like streamer patterns. 
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Fractal properties are common to the dielectric breakdown phenomena H] which range from an atmo- 
spheric lightning [|],[|J. [|| to electric treeing in polymers |5] and are of scientific and technical importance 

Although the actual physical processes can be quite different in these phenomena, the global prop- 
erties of the resulting discharge patterns are very similar. Filamentary gas discharges on insulating surface 
exhibit remarkable similarities to breakdown phenomena in long gaps e.g., to atmospheric lightning , and 
thus offer the possibility to perform well-defined model experiments in laboratory. The application of sophis- 
ticated diagnostic tools, such as streak cameras, high-speed oscilloscopes, and time-resolved spectroscopy ||, 
in surface discharge experiments has improved the basic understanding of how a highly conducting phase, 
the filamentary "leader" channel, advances into a nonconducting medium as the surrounding gas. 

The surface discharge in compressed SF e gas have been studied in detail by Niemeyer and Pinnekamp|i"o|. 
The parameters were controlled in such a way that the experiment produces, to a good approximation, an 
equipotential channel system growing in a plane with a radial electrode from a central point. The experiment 
shows that the dielectric breakdown pattern has a fractal structure. 

The stochastic model containing the essential features of the fractal properties of the dielectric break- 
down was introduced by Niemeyer, Pietroniero, and Wiesmann (NPW) The introduced stochastic 
model made it possible to simulate the growth of a fractal structure which resembles a Lichtenberg figure Q . 
The basic assumption of the dielectric breakdown model (DBM) was that the growth probability depends on 
the local field. This model naturally leads to fractal structures. In |ll| the Haussdorf dimension and other 
fractal properties of the structure are determined and their close relations to the other fractal structures e.g., 
DLA (the diffusion limited aggregation) [Q, is shown . 

The DBM method was generalized by Wiesmann and Zeller [jl3| in two ways. Firstly they introduced 
a critical field for growth F c . The growth probability was assumed to be proportional to the local field Fi oc 
if Fioc > F c and zero if Fi oc < F c . Secondly they have introduced an internal field in the structure F s . The 
potential in the structure was no longer equal to the potential Vq of the connecting electrode but equal to 
Vq + F s ■ s, where s is the length of the path along the structure which connects the point to the electrode. 

In the recent experiment |3) a transverse high magnetic field was applied during the discharge 
evolution and thus any spatial restriction of the surface discharge was avoided in order to use a locally 
sensitive probe. 

When a rectangular high-voltage pulse of about 500ns length, 20kV amplitude, and less than 10ns 
risetime is applied to a point-to- plane electrode system, the electrodes being separated by a thin dielectric 
film ( thickness i=s lOOum) , a discharge propagates in the gas just above the surface of the film. The discharge 
pattern was recorded directly with high resolution by using a photographic film as dielectric plate. 

The experiment shows the spatial evolution of a negative surface discharge in a nitrogen atmosphere as 
a function of the magnetic field ( Fig. 2 of Ref. ||). At B = a very bright starlike pattern develops. At 
moderate magnetic field (up to IT ) , the leader channels are bent and appear to have a circular shape outside 
the central electrode region. 

The radius of curvature is of the order of 1cm at IT. The direction of the bending corresponds to the 
movement of electrons in crossed electric and magnetic fields. With the increase of the magnetic field 
the radius of curvature decreases, the channels approach each other and branching sets in. At the highest 
applied magnetic field of 12T, circular-shaped current filaments are only found in the outmost regions of the 
discharge pattern where they can develop undistributed by the fields of neighboring leader channels. 

One can summarize the experimental results as follows: 

(1) The observed bending effect cannot be related to the movement of a single charged particles in crossed 
electric and magnetic fields, which should result in the curvatures of the order of magnitude of 15/im at IT 
for the electrons, i.e., the Larmor radius r c — v c /uj c of the cyclotron orbit of an electron withe enough kinetic 
energy to ionize gas molecules. It has be assumed that the macroscopic value of the radius of curvature results 
from the drift movement of the Larmor centers of the gyrating electrons. The average radius of curvature 
follows a power law, R oc B~ a . The exponent a depends on the nature of gas. 

(2) The circular shape of the current filaments is not expected to occur for our point-to-plane electrode 
geometry. 

(3) The complexity of the discharge pattern can be controlled by the magnetic field strength. It appears 
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that an increase in complexity is induced on the scale of the bending observed, and not on the microscopic 
scale of the cyclotron radius. 

In what follows we modify the model introduced by Wiesmann and Zeller Jl3j for the case of external 
magnetic field. 

Let us consider a two-dimensional square lattice in which a central point represents one of the electrodes 
while the other electrode is modeled as a circle at large distance, which represents the geometry of the 
experiment. The discharge starts at the central electrode and grows by one lattice bond per growth step. A 
bond connects a lattice point with one of eight adjacent lattice points as it is described in Fig.l. Once a given 
point is connected to the discharge structure by a bond, it becomes part of the structure. The potential of 
the central electrode is 0, the potential of a point in the structure is Vi t k = Vi >rn + Vr * I, where Vi >m is the 
potential of a points from which growth go on, Vr is an internal field in the structure and I is 1 for bonds 
parallel to the grid and is \/2 for diagonal bonds. 

The growth is computed as follows: First the Laplace equation is solved with the boundary conditions 
determined by the electrodes and the discharge structure. Then the local field Fi oc between a point which 
is already a part of the structure [i, fc) and a new adjacent point («', k') is calculated: 



Fjocft M', *') = (!) 

where tp t l is the solution of Laplace equation. The breakdown can occur only if the local field is greater 
that the critical field of growth Fc- The probability that a new bond will form between a point which is 
already a part of the structure and a new adjacent point p(i, k — > i' , k') is calculated as a function of the 
local field Fi oc between the two points: 



F loc {i,k,i',k') < F c 

(2) 

Fi oc (i,k,i',k') > F c 

where a power-low dependence with exponent r\ is assumed to describe adequately the relation between the 
local field and the probability. The sum in the denominator refers to all possible growth processes. A new 
bond is chosen randomly with probability distribution (|^) and added to the discharge pattern. With this 
new discharge pattern one starts again. More detailed description of this model is done in jl3| . 

Let us take into account the magnetic field. The magnetic field changes the probability distribution (||). A 
moving particle in the magnetic field experiences the Lorentz force Fl ~ V* x H which acts perpendicularly 
to its velocity. Consider each step of growth like a superposition of two processes. The first step is choosing 
a new bond using the probability distribution (||) , and the second step is taking into account the probability 
of deviation of the growth due to the magnetic field, pn- 

If after first step of growth, for example, the bond from the dot to the dot 4 was chousen (Fig.l), the 
growth will occurs to the dot 3 with the probability pu and to the dot 4 with the probability 1 — pn ■ The 
new probability of growth can be written as 



p(i, k — > i', fc') = 



., = p(i, fc -> i', k'){l - PH (t, fc, i', fc')) + pis k -> i", k") PH (i, k, i", k") 
P[h ' T,(p(hk^i',k')(l-p H (i,k,i',k'))+p(i,k^i",k")p H (i,k,i",k")) U 

where the point (i", fc") is the neighboring point with respect to the point {if , fc') in the clockwise direction 
with respect to the point (i, fc) and the sum in the denominator refers to all possible growth processes. The 
probability [| was used in computer simulations in place of the probability || with the same algorithm. 

Let us turn to the probability pu of deviation of the growth due to the magnetic field. In our model, 
during the process of growth, two constant forces act on charge cariers, Fl and Fh (Fig.l). When the 
resulting force is near the dot 3, pn — > 1; when the resulting force is near the dot 4, pa — > 0. It is clear 
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that the probability pu is proportional to ■^ L - . If Fi oc (i, fc, i' , fc') < Fq then the growth does not occur and 

Ph should be zero. If -jp^- > 1 then we let pu = 1- The Lorentz force Fl is proportional to the velocity of 
charge carrier. In our model we can take into account only a local velocity, which arrises in the first step of 
growth process due to the acceleration in the local field Fi oc . So, the Lorentz force should be proportional 
to y/lFi oc (i,k,i' \k'). 

Based on these considerations, we choose the probability of deviation of the growth due to the magnetic 
field, ph, in the following form: 



F loc (i,k,i',k') < F c 



PH(i,k,i',k') = < 



F L (i,k,i' ,k') 
Fioc (i,k,i' ,k') 



-g$^<i, F lo s^',k')> Fc 



(4) 



1 . > l > Fioc(i,k,i',k')>F c 

where Fl(i, k, i', fc') = y/lFi oc (i, fc, V , k')H and H is the value of the magnetic field. 

In our computer simulations we consider a 500 x 500 lattice. The solutions of the Laplace equation were 
obtained by the iteration method Before starting of each realization of growth we performed 20000 

iterations and after each step of growth the number of iterations was 40. This procedure gives a good 
convergence. The number of particles in clusters was 9000. 

The fractal dimension was calculated by the method described in [fll|| . For every realization we plotted 
the log N (R) versus log R where N (R) is the number of particles belonging to the structure and being within 
a circle of radius R. The fractal dimension is obtained by fitting a straight line to the data scaling region. 
For every set of the same parameters of the model (H, Fc, Vr, rf) we made about 100 realizations. Thus 
the statistical fluctuations were reduced. 

Let us discuss the results of these calculations. 

We start our simulations with the case of the zero magnetic field and the zero values of the parameters 
Fc and Vr in order to compare our results with the results of the different authors. In the Table. 1 we present 
the dependence of the Hausdorff dimension D on the exponent parameter r\. 



V 


D (our results) 


D (according to |fll|) 


D (according to ]13| ) 


0.5 


1.89 ±0.02 


1.89 ±0.01 


1.92 


1 


1.73 ±0.02 


1.75 ±0.02 


1.70 


2 


1.6 ±0.03 


- 1.6 


1.43 



Table 1. Dependence of the Hausdorff dimension D on the exponent r\ used in the relation between 
probability and local field (Eq.||). 

Our results are in a good agreement with the results, obtained earlier. 

The example of the computer-generated discharge pattern (Lichtenberg figure 0) corresponding to the 
following set of parameters: H = 0, i) = 1, F c = and Vr = is shown in the Fig. 2. In the Fig. 3 we show 
the computer-generated discharge pattern in the presence of the external magnetic field. The white lines 
in the figures 1 and 2 correspond to the leader channels. Unlike the pattern, presented in the Fig. 2, the 
leader channels in the magnetic field are distorted and appeare to have a circular shape outside the central 
electrode region. The direction of the bending corresponds to the movement of electrons in crossed electric 
and magnetic field. At the lower left corner of the Fig. 3 we plotted the cyclotron orbit of an electron. One 
can note that the Larmor radius of an electron is about two orders of magnitude smaller then the radius of 
curvature of the leader channel. 

The saturation of the fractal dimensionality, with growing magnetic fields, at the value of D = 1.67 , 
which is one of the main results of this paper, is presented in Fig. 4. The plot starts from the value of 
D = 1.65 , in the absence of magnetic field, which is smaller than the fractal dimensionality reported in JlT[] . 
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Such difference results from the fact that in the critical field value for the breakdown was not taken into 
account. We have improved their calculations by introducing the minimal value of the electric field for the 
breakdown between two successive points. In this case the breakdown pattern is more directionally selected, 
and a lower fractal dimensionality results. With the increase of the external magnetic field, the MFD growth 
and finally saturates at an universal for high magnetic fields values of 1.67 ± 0.01. The growth of MFD with 
the magnetic field could be expected, since the curved trajectories fill up the space more densely than the 
straight ones. The existence of an universal limit, hovewer, is far from being obvious. Following the directed 
percolation models, one could think that the saturation of MFD will occur at D = 2 . Fig 4. shows clearly 
that in this system the MFD saturates due to the physics of current carrying streamers. 

An unexplained feature in the experiments of || is the 'spider-legs' form of the breakdown pattern in the 
absence of the external magnetic field. We outline here that the streamer currents are rather strong, 104- 100 
A, and their influence on the streamer pattern can be very important. To describe this phenomenon we have 
taken into account the magnetic interaction between current carrying streamers, in the framework of the 
modified active walker model (MAWM) which is described in what follows. The results of our calculations, 
presented in Fig. 5, show that our model correctly describes the experimentally observed 'spider-legs' effect. 

Active walker models have been used to describe different pattern formation problems |lj, |1| . In these 
models the walkers moovement is subject to the influences of the environment and vice versa. We describe 
the leader channel propagation in the magnetic field using the active walker model. The Lorentz force acting 
on the fast-moving electrons is particularly effective in the high-field regions in the leader tips, where the 
channel formation takes place. 

Let us consider again a two-dimensional square lattice in which a central point represents one of the 
electrodes while the other electrode is modeled as a circle at large distance. The discharge starts at the 
central electrode, so initially several walkers are set in the vicinity of it. The walkers move in a potential 
which is the solution of the Laplace equation with the boundary conditions determined by electrodes and 
discharge structure. During a step of growth each walker moves. The solution of Laplace equation is found 
by iteration method |ll| after each step of growth. 

When a walker moves to a point, this point starts belonging to the breakdown structure. The potential 
of the central electrode is 0, the potential of a point in the structure is Vi t k — Vi >m + Vr * I, where V% m is the 
potential of a points from which growth go on, Vr is an internal field in the structure and I is 1 for bonds 
parallel to the grid and is v2 for diagonal bonds. In the absence of the magnetic field, the probability of a 
walker step is a function of the local field F\ oc |l4| . The breakdown occurs only if the local field is greater 
then the critical field of growth Fc- 

The magnetic field is taken into account by the following way. We add to the local field in the direction 
perpendicular to the previous move of the walker the Lorentz force which is proportional to the magnetic 
field. The magnetic field acting on the i walker is Hi = Hq + H}, where Hq is external magnetic field and 
Hj is the field created by currents of the breakdown structure. We calculate H\ by means of the Biot-Savart 
law: 



where / is the current in a channel, d~l?k,i is I element of the k channel and ~r* is the vector pointed from 
the position of d~s*k,i to the position of the i walker. 

The main results of this part of studies are presented in Fig. 5. 

To summarize, we have generalized and modified the existing dielectric breakdown models to explain the 
experimental observations of the propagation of a streamer near an insulating surface under the influence 
of a transverse magnetic field. We have introduced the concept of the Magnetic Fractal Dimensionality 
(MFD) and have obtained its saturation with growing magnetic fields. The Universal Magnetic Fractal 
Dimensionality (UMFD) equals 1.670 which is superior to 1.65, the one in the absence of a magnetic field. 
Inclusion of the magnetic interaction between the current-carrying streamers results in the 'spider-legs' like 
streamer patterns at lower fields, which corresponds to the experimental observations. 




(5) 
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Figure Captions 



Fig.l. 

Bonds connecting a lattice point with one of eight adjacent lattice points in a two-dimensional square 
lattice. The geometry of the experiment is represented by a central electrode at some point while the other 
electrode is a circle at large distance. The discharge starts at the central electrode and grows by one lattice 
bond per growth step. 

Fig.2. 

Computer-generated discharge pattern (Lichtcnbcrg figure) corresponding to the following set of param- 
eters: H = 0,t] = 1,F C = and Vr = as explained in the text. The white lines correspond to the leader 
channels. 



Fig.2. 

Computer-generated discharge pattern in the presence of the external magnetic field. The leader channels 
in the magnetic field are distorted and appeare to have a circular shape outside the central electrode region, 
due to the action of the magnetic field. 

Fig.4. 

The saturation of the fractal dimensionality, with growing magnetic fields, at the value of D = 1.67 . 
Fig.5. 

The 'spider-legs' form of the breakdown pattern folowing from the magnetic interaction between the 
streamer currents: a) current-current interactions arc not took into account, H = 0; b)current-current 
interactions are took into account, H = 0; c) H ^ 0. 
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